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Abstract 

The Chapman-Enskog method of solution of the relativistic Boltzmann equation is generalized 
in order to admit a time-derivative term associated to a thermodynamic force in its first order 
solution. Both existence and uniqueness of such a solution are proved based on the standard 
theory of integral equations. The mathematical implications of the generalization here introduced 
are thoroughly discussed regarding the nature of heat as chaotic energy transfer in the context of 
relativity theory. 
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I. INTRODUCTION 



The standard treatment of the Boltzmann equation (BE) provides microscopic support 
to the transport equations commonly used in science and engineering problems QQl. It 
is well known that whereas the local equilibrium solution to the homogeneous BE leads to 
the Euler equations of hydrodynamics, its first order correction in the Knudsen parameter 
yields the Navier-Stokes-Fourier transport equations of fluid mechanics for a dilute gas. In 
the standard non-relativistic method of establishing the first order correction to the local 
equilibrium function, all partial time derivatives are replaced with spatial gradients using 
the Euler equations. Both existence and uniqueness of the solutions built in this scheme 
have been shown using the theory of integral equations 

On the other hand, in relativity theory time is simply a new coordinate ct = x 4 . According 
to this fact, the direct application of the conventional method used to obtain the first order 
correction to the local equilibrium function is shown here to lead to a first order correction to 
the distribution function in terms of both spatial and time derivatives of the local variables 
namely, the generalized four-component thermodynamic forces. This is done following the 
tenets of relativistic linear irreversible thermodynamics as shown in Ref. |4h We wish 
to emphasize that in this work the Meixner version of non-equilibrium thermodynamics is 
developed. However, in most standard works on relativistic kinetic theory [f], Q, [|| these 
ideas are ignored. The consequence of the absence of time components in the structure of 
the thermodynamic forces is on the one hand, the obtention of transport equations which 
apparently violate causality, and on the other hand the tensorial structure of the heat flux in 
relativity, which is still a subject of debate [9||. The main motivation of this work is precisely 
to critically analyze these fine points. 

The mathematical procedure here used to show that the solution obtained, considering 
four component thermodynamic forces, exists and is unique is based on the standard theory 
of integral equations. One requires orthogonality between the inhomogeneous part of the 
equation and the collisional invariants. This relation implies the validity of Euler equations. 
Nevertheless the use of these equations is not required for the construction of the solution 
to first order in the gradients. These results have strong implications on the structure of 
the four component heat flux as will be discussed in this work. Finally, the usual subsidiary 
conditions are invoked in order to adjust integration constants. This has, as will be discussed 
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in this work, strong implications on the structure of the heat flux tensor. 

The paper is divided as follows. In Section 2 the relativistic BE and some elements of 
the Chapman-Enskog method of solution are briefly reviewed. The general form of the first 
order in the gradients solution to the relativistic BE is proposed in Section 3 and both its 
existence and uniqueness are proved in Section 4. Finally, Section 5 includes a discussion of 
the implications of the results and some concluding remarks. 

II. THE CHAPMAN-ENSKOG EXPANSION 



The starting point of the generalized formalism is the special relativistic Boltzmann 
equation for a simple system in the absence of external forces |5|, Ifjj: 



v a f, a = j(f,n , a) 

which describes the evolution of the single particle distribution function / = 
f(x u , v u \ n(x u , t), u(x u ,t), e(x u ,t)) namely, the molecular number density in phase space. 
In Eq. (Dp) 



df_ 

dx k 

c dt 



(2) 



and the index a, as well as all greek indices in the rest of this work, runs from 1 to 4 while 
the latin ones run up to 3. The local variables, number density n, hydrodynamic velocity 
u k and internal energy e per unit of mass are defined in terms of / as follows 

// / fodv* , (3) 



nu 



f^v k dv* 



ne = mc / fjv dv* . 



(4) 



(5) 



Here v M is the molecular velocity four- vector which we denote in terms of the three velocity 
w as 

7c 

where 7 = (1 — w 2 /c 2 ) is the usual relativistic factor. Also, the differential velocity 
element can be written as {3| 



(6) 
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Cu ID ) 

dv* = 7 5 --f- = Aire 3 y/i*^ld 1 , (7) 



and, as usual, c is the speed of light. Equation (pQ) implies that changes in the distribution 
function are due to collisions, represented on the right side of Boltzmann's equation through 
the so called collision kernel J (/, /') which is given by 

J (f, f)= f I (/' fi-f h) * gcMdvt . (8) 



In Eq. (jHJ) primes denote quantities after a binary collision between particles with velocity 
v and V\ takes place; o (Q) dQ is the differential cross section element and g is the relative 
velocity. 

To solve Eq. (pQ) following what is now known in the literature as the Chapman-Enskog 
method jl], originally due to Hilbert, the distribution function is expanded in a power series 
of the Knudsen's parameter e around the local equilibrium distribution function 



f = /W (1 + e0« + eV 2 ) + -..). (9) 

The expansion parameter e is a measure of the relative magnitude of the gradients of the 
local variables within a mean free path and the characteristic size of the system. The local 
equilibrium distribution function, solution to the homogeneous relativistic BE J (/, /') = 0, 
is the well known Jiittner function. For particles of rest mass m, relativistic parameter 
z = ^2 and in the non-degenerate case [si, o|: 

As usual, k is the Boltzmann constant, n and T are the local density and temperature, 
respectively, and /C 2 represents the modified Bessel function of the second kind. 

Substitution of Eq. ((9)) in Eq. (pQ) and considering only linear deviations in e from local 
equilibrium yields 

v a f iO ) = j( / (O) (1 + e0))/ (O) (1 + e0O ) j (U) 

where ecj)^ = <p and f^°\ as pointed above, is the solution to the homogeneous equation 
namely, 

J(/ (0) ,/ (0) )=o, 

or 

f (o) f m = f{0)f (o) _ (12) 



Introducing the relation given by Eq. (112!) in the right hand side of Eq. (TTT1) leads to 

J (/, /') = / {0) f [ Wi + <f>' - h - <P] fi% (O) ^cKMt^ . (13) 



Thus, this procedure yields a linear inhomogeneous integral equation for 0, namely 

= /0)C (0) , (14) 

where 

C (0) = / / {0'i + 0' - 0i - 0} (SI) gdtldvl , (15) 



is the linearized Boltzmann collision kernel. The following section is devoted to the obtention 
of the solution to Eq. lfl~4|) . 



III. SOLUTION TO FIRST ORDER IN THE GRADIENTS 

Following the theory of integral equations as presented in Ref. js]], the general solution to 
Eq. ( fT4l) may be written as the sum of a solution to the inhomogeneous equation 0p and 
an arbitrary linear combination of the solutions of the homogeneous equation 0#: 

= 0iT + 0P • (16) 

The homogeneous equation has five solutions, the collisional invariants namely, 

1=0. (17) 
rrvyv I 

Notice that, for \i = 4 Eq. (fT7l ) corresponds to the particle number conservation. 

A solution given by Eq. ( 1T61) exists if these invariants in 0# are orthogonal to the left hand 

side of Eq. ([Till namely 

mv 
mjv 

The relations in Eq. ( fTSl correspond precisely to the relativistic Euler equations and thus 
the orthogonality conditions are satisfied. However we recall that there is an infinite number 
of solutions since 0# is an arbitrary linear combination of collisional invariants. The unicity 
is achieved when the arbitrary constants are determined with the use of the subsidiary 
conditions as we shall argue below. 

In order to find a particular solution to the inhomogeneous equation, we start by calcu- 
lating the left hand side of Eq. ( fT4l) . In the comoving frame and assuming that no external 
forces are present 



,4 



v a f ( Vdv* = . (18) 
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, Q/ (0) 



(19) 



V <9n ' <9T 

Since a = 1, 2, 3, 4 each gradient on the right side of Eq. ( fT9l ) has four components. It is 
precisely at this point where our formalism deviates from the standard one 1, 5, f|. In con- 
ventional relativistic treatments, one separates time and space derivatives using a projection 
operator and introduces Euler's equations in order to eliminate the time derivatives of the 
local variables. However, special relativity is a four dimensional formalism in which time 
corresponds to a coordinate. Because of this feature we here propose that the four-gradients 

n 

are to be considered with the same status as the generalized thermodynamic forces [10]. 
This is why we do not resort to the Euler equations to eliminate time derivatives. Also, 
in most standard relativistic treatments, the temperature gradient is written in terms of 
the pressure gradient in order to introduce, through the momentum balance equation, an 
acceleration term so that an Eckart-type constitutive equation is sustained [5j. However, in 
the Meixner-type scheme here considered, such a substitution needs not to be carried out 
since the thermodynamic forces are T v and n „ exclusively. 

In order to proceed, we first notice that since the calculations will be performed in the 
comoving frame, = and = 0. Thus, using the equation of state for an ideal gas one 



can write for the space components of the gradients that 

n T 

For the time component of the density gradient, in the comoving frame, one has 

n 4 = , 

since the continuity equation is satisfied to any order in e. Introducing Eqs. ( 1201 ) and ( l2Ti l 
in Eq. (TT9l) and calculating explicitly the derivatives of the Jiittner function, the integral 



(20) 



(21) 



equation for (j)p can be written as 



v k Ti 



zT 



* f(0) 



7 - Q 



v 4 T, 



T 



£ /( o) 



7 



/ (0) C (<f> P ) 



(22) 



where Q (I) = /C 3 (±) /JC 2 (i). 

Since the thermodynamic forces present on the left hand side of Eq. (l22l) are independent, 



we seek a solution as a linear combination of such forces: 

0p = , 



(23) 



which is a direct generalization of the structure proposed in the non-relativistic case. The 
coefficients A" ', which are in general functions of v a , n and T, are solutions to 
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IT 
Z 



l - Q 



C (^*) , 



for the first three components and 



- 7-0 - 



+ 1 



/<°>c (^ 4 ) 



(24) 



(25) 



for the fourth component. It now remains to show that one may obtain a unique solution 
to Eq. jHD by using the subsidiary conditions. This task will be performed in the following 
section. 



IV. SUBSIDIARY CONDITIONS AND UNIQUENESS OF THE SOLUTION 



According to the previous discussion, the general solution to the relativistic BE, Eq. (DQ), 



to first order in the Knudsen parameter, has the general form 

T 

(p = A u ^r + a + a v v v , 



(26) 



where the first term corresponds to the particular solution, Eq. ( 1231 ). while the remain- 
ing ones, as mentioned above, are an arbitrary linear combination of the solutions to the 
homogeneous equation, namely the collisional invariants. 

The components of A u are given by the solution to Eqs. (1241) and (l25lh However, in order 
to uniquely determine the complete solution 0, the constants a and must be determined. 
This is usually achieved by arbitrarily assuming that the local macroscopic variables are 
determined only from the local equilibrium state. That is, 



r(0) 



rwyv 



dv* = 



(27) 



A" 



(28) 



which are the so-called subsidiary conditions [l], Q). Since spatial and temporal terms in <pp 
have opposite parities (in v), it is convenient at this point to explicitly separate such terms 
by proposing that 

a' (7) v k 
a" (7) v 4 

Thus, Eq. ( |27l ) can be written as 

J (a' (7) v k T k + a" (7) v A T A + a + a v v v ) ijjf^dv* = , (29) 

where if) are the collisional invariants and 7f 4 . For if) = v 4 , the odd parity in a' (7) v k T^ 
and ctkV k yields zero values for the respective integrals and thus 
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J (a" (7) v 4 T A + a + a 4 v 4 ) 7 / (0 W = . (30) 

In a similar way, for ift = v 1 the only non-vanishing terms due to parity are a' (7) v k T tk and 
a k v k , whence 

J (a' (7) v k T k + a k v k ) v e f^dv* = , (31) 

which is non-zero only for £ = k. In this case 

J {a' (7) v k T k + a k v k ) v £ f^dv* = ~ J {a 1 (7) T k + a k ) ^w'f^dv* (32) 

Thus, the three subsidiary conditions for if) = v e may be written as 

j (a' (7) T k + a k ) -fvffWdv* = . (33) 

while the subsidiary condition corresponding to the fifth collisional invariant, if) = 4 , is 

J {a" (7) v 4 T 4 + a + a 4 v 4 ) 1V 4 f {0) dv* = . (34) 

From Eq. f|33l) . it is straightforward to obtain the relation 

^ ~ ~ T > k J>yWf«»dv* ' (35) 

which implies a k oc T k . Thus, one can redefine the coefficients appearing in Eq. (1291) as 
follows 

a' (7) v k T k + a k v k — a' (7) u fc T fc , (36) 

This implies that the subsidiary condition given in Eq. (13T1) now reads 

a' (7) 7V/ <°W = . (37) 



On the other hand, from Eqs. ( l30l) and ( l34l ) an inhomogeneous system of equations for a 
and 04 is obtained: 

g n a + gi 2 a 4 = s{T A , (38) 
9*\ol + 922 ®4 = s 2 T A , (39) 

where the coefficients are 

jfWdo* g 12 = g 2 i = I j 2 cf^dv* g 22 = [ cWf {0) dv* , (40) 



0u = 

and the coefficients in the inhomogeneous terms are given by 



si 



a" (7) -fcf^dv* s 2 = - a" (7) c 2 7 3 / {0 W , (41) 



The solution to this system of equations is 



^ sig 22 - s 2 g 12 
a = T 4 , (42) 

gng22 — gug2i 



rp s 2 gu - giggi f 
a 4 = -/ i 4 . (43) 

9ll922 — 912921 

Since, from Eqs. (1421 ) and ( 1431 ). both a and a 4 are proportional to T 4 one can redefine 

a" (7) t» 4 T 4 + a + a 4 t> 4 — ► a" (7) t» 4 T 4 , (44) 

and therefore the two corresponding subsidiary conditions are now given by 

a" (7) l 2 f {0) dv* = , (45) 



a" (7) 7 3 / ( °W = , (46) 

Putting these results together, the unique solution to the relativistic BE to first order in 
the gradients according to the Chapman-Enskog expansion has the general structure 

/ = / (0) (l + ^) , (47) 

where A u is defined in Eq. (1251) and the coefficients a' (7) and a" (7) are given by the solution 
of the integral equations 

v k 



z 

4 f (0) 



/ (0) (7-$Q))=C(a'( 7 )t; fc ) 



(48) 



C (a" (7) v A ) , (49) 



subject to the subsidiary conditions given in Eqs. ( 1371 ). (1451) . and ( 1461) . 
Notice that, in the non-relativistic case as well as in the relativistic treatments given in the 
literature g Q, the absence of a f o Ur th co mp o„e„t of ^ !ea ds ,o a ao m o g e„eo US SyS te m 
instead of Eqs. (|38l) and (|39l) and is thus responsible for the appearance of the trivial 
solutions a = <S 4 = 0. The consequences of this substantial difference are discussed in the 
next section. 



V. DISCUSSION AND CONCLUDING REMARKS 

By considering the four-dimensional gradients of local variables as thermodynamic forces 
within the Chapman-Enskog method of solution to the BE, a generalized first correction to 
the equilibrium distribution function has been found in the context of Meixner's relativistic 
irreversible thermodynamics. In the comoving frame, this procedure leads to a first order 
in e solution with the general structure given by Eqs. (1471) and (f28|) . The proportionality 
coefficients, functions of the local variables and 7, can be obtained by an appropriate decom- 
position in orthogonal functions which simplifies the integral equations (1451 ) and (1491) . Thus, 
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the solution here obtained is suitable for the calculation of relativistic transport coefficients 
consistent with the phenomenological formalism proposed in Ref. Such a procedure is 
beyond the scope of this work and will be addressed elsewhere. 

However one important point should be discussed here, namely the definition of heat 
flux in the context of relativity. We recall that in the non-relativistic theory this quantity 
is related to the transport of chaotic kinetic energy in agreement with the interpretation 
proposed by Clausius in 1857 fll|. The question here is: can this definition be extended 
within a relativistic framework? This, in principle, may be achieved by resorting to the 
relativistic version of the Maxwell-Enskog transport equation. Its explicit form was derived 



seven years ago 



121 ] by two of us. This equation, in the notation here introduced reads 



J,% = «' / v"fjdv- 



n{^-))^0 (50) 

where ip = mv M , m^v A is a collisional invariant and the bracket represents the statistical 
average defined, for any function g (7), as 

U) = \J f9ldv* (51) 

For g = mv 4 the heat transport equation is obtained. In particular, the heat flux will be 
given by 

'[Q] 

This is a rather striking result because the fourth component of this vector turns out to 
be precisely the energy density as defined in Eq. (JSj) - This implies that, because of the 
subsidiary conditions, such component will be identically equal to zero for any perturbation 
0M with n > 1. Thus we recover the usual form which reduces to the equation obtained by 
other authors resorting to the use of a spatial projector If], Q, 

Indeed to clarify this statement, we recall that Eckart proposed the following stress-energy 
tensor (§} for a simple, non-viscous fluid 

TIE I I 

T^ u = -^u^Uy + php V + -^u^Jy^y + -^u u J [Q] ^ (52) 

where h^ v = g^ v — u fJ- u u /c 2 is the projector introduced by Eckart in order to eliminate 
dissipation in the time-axis. This work thus provides further support for this assumption. 
However, T^ u as in Eq. ( 1521) is in principle coupled to the Einstein tensor by means of 
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the field equation: 

= KT^ (53) 

where K, is the coupling constant. The point here is how heat would affect the gravitational 
field if it is placed in the stress energy tensor. Indeed, heat is identified with molecular motion 
and its capability of transforming itself to mechanical energy by any process is restricted. 
On the other hand, the introduction of J[q\ v in Eq. (l53l) seems to allow a transformation of 
heat to work, via the gravitational field, with unitary efficiency. To the authors' knowledge, 
direct attempts to obtain solutions to the system (l53|) in the presence of heat and to analyze 
the dynamical consequences of the solutions are practically absent in the literature. 

An alternative way of introducing heat in relativistic irreversible thermodynamics consists 
in a generalization of Meixner's approach |4|. In it, heat is not included in the stress 
energy tensor but introduced directly in the total energy flux. The construction of the 
entropy production, consistent with Clausius' idea of uncompensated heat, and consequent 
enforcement of the second law of thermodynamics in this context leads to the proposal of 
the constitutive equation 

Jfa = -kT* (54) 

where k is the thermal conductivity. Equation ( 1541) coincides with the constitutive equation 
in Eckart's formalism in the comoving frame where the acceleration term vanishes. 

It is important to emphasize at this point that for the purpose of this discussion, the 
key difference between both formalism is the role of the heat flux. On the one hand it is 
included in the stress energy tensor in Eckart's formalism, as pointed out above. On the 
other hand, it is conceived only as energy in transit in the Meixner-like formalism and is 
introduced separately in the total energy flux. Clearly the latter conception does not predict 
any effects on the geometry of space-time due to heat since is not included in and 
thus the Einstein field equation remains unchanged. An experiment that could clarify this 
difference has already been suggested in Ref. 

On the other hand, the validity of Eq. (I54ll should be sustained by solving the relativistic 
Boltzmann equation. The calculations performed in the previous sections show that this is 
not the case. Indeed, the heat flux has no fourth component which is in agreement with 
other theories but for entirely different reasons. 

We believe that this poses a rather intriguing question. Is it possible or not to think 
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of heat as a chaotic molecular motion in a relativistic framework? The correct answer to 
this question is crucial when considering both the structure of transport equations and the 
consistency of any relativistic kinetic theory with the second law of thermodynamics J|. 
These questions merit further study. 
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